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COHERENT COMPOSITION OF SIGNALS BY MEANS OF PROGRESSIVE 
FOCAL CORRECTION 

TECHNICAL FIELD 
5 Physical Technology 

Electronic system for clinical echographies and non-destructive testing 

Other applications: radar, sonar, seismology. * 

PRIOR ART 

10 In various fields of application, a set of N transducers or array, receives signals 
from one or multiple sources located at a certain distance and orientation from 
the array. This is a common situation in radio astronomy, radar, sonar, 
seismology and, particularly in the field of array-based ultrasonic technology 
(clinical or industrial echographies, in this case as applied to Non-destructive 

15 testing), where N is, typically, in the range of 32 to 256 or more. 

In these situations, each of the N elements of the array receives the signal 
generated by each source with a certain temporary delay in relation to the other 
elements. The function of a beam former is to obtain one sole signal from all the 

20 N signals received in which the temporary delays for each possible distance at 
a given orientation are compensated, increasing the discrimination capacity of 
the system (greater resolution, signal/noise ratio and contrast). By repeating this 
process in a variety of directions by means of linear or angular scans, sufficient 
information is obtained to compose a 2D or 3D image from the source signals. 

25 In the specific case of ultrasounds, the N signals received are echoes produced 
by the discontinuities in the media that have been sonified by an emitter that 
may be the same receptor array (or a subset of same) or other transducer 
(array or monoelement). 

30 In the reception end, it is desirable to compensate the delays at all the points or 
foci along the propagation line in a given direction. This requires the continuous 
modification of the time delays to be able to adapt them for the incoming echoes 
produced at dynamically variable distances from each element. This dynamic 
focusing technique provides the greatest resolution, contrast and signal/noise 

35 ratio. 
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Since it is technically impossible to modify delays in a continuous fashion it is 
essential that they are discretized with an excellent temporal resolution for each 
array element and for each distance, typically in the order of 1/16 to 1/64 of the 
5 fundamental signal period. 

This fact has spurred the search for dynamic focusing reception methods that 
are simple and efficient and that have been adapted to the available technology 
through time. In this manner, techniques based in CCD devices actuated by 

10 clocks that determine the delay time by their period have been described 
[Macowski, "Ultrasonic electronic lens with reduced delay range", US Patent N° 
4058003, Nov, 1977], of in switchable delay lines [Maslak, , US Patent N° 
4550607, "Phased array acoustic imaging system" Nov, 1985]. These methods 
introduce switching artifacts, are bulky and provide limited resolution, and are 

15 practically in disuse. 

Another alternative, also analogic, is based in the use of mixers and phase 
shifters that effect the temporal delays in the frequency domain [Malask et al., 
"Dynamic focused linear phased array acoustic imaging system", US. Pat 
20 4699009, Oct, 1987]. One the problems associated to this type of techniques is 
how difficult it is to adapt them to transducers operating at different frequencies. 

On the other hand, the use of phase scramblers based in the CORDIC 
algorithm [ R. Andraka, "A survey of CORDIC algorithms for FPGAs", Proc. 

25 1998 ACM/SIGDA VI Int'l symposium of FPGAs, Monterrey, CA, pp. 191-200, 
Feb. 1998], and on their application to dynamic focusing can be found in [ 
O'Donnel et al., "Method and apparatus for digital phased array imaging", US 
Patent N° 4983970, Jan. 1991] and in [Engles et al. "Dynamic phase focus for 
coherent imaging beam formation", US Patent 511695, May 1992], in which 

30 there is significant hardware involved. 

.» 

The availability of the VLSI technology allowed the formulation of new proposals 
based in delaying the signals once digitalized. Each element of the array is 
associated to an independent acquisition and processing channel (focusing). 
35 The N channels outputs are added to obtain the resulting signal (coherent 
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summation). To operate in real time, each channel has two memory spaces: 
one for data, where the acquired samples are stored, and the other assigned to 
focusing to establish the delays to be applied. 

5 In principle, delay resolution was limited to the signal sampling period itself, as 
established by the Nyquist criterion, between approximately 1/3 to 1/5 of the 
signal's period, using the variable length shift registers [Welles II et al., 
"Architecture for ultrasonic imaging", US Patent 4796236, Jan 1989]. The low 
temporal resolution produces delay quantification lobules that reduce the 
10 dynamic range, contrast and the signal/noise ratio of the image [D.K. Peterson, 
G. s. Kino, "Real-time Digital Image Reconstruction: A Description of Imaging 
Hardware and an Analysis of Quantization Errors", IEEE Trans. Sonics 
Ultrasonics, 31, pp. 337-351, Jul 1984]. 

15 The option of using a greater sampling frequency is not advisable due to the 
greater cost of the components, energetic consumption and the volume of data 
to store and process. To solve this problem several techniques have been 
proposed that can be grouped in two kinds: 

20 1. Gross delay with a temporal resolution equivalent to the sampling period, 
and a fine delay achieved by interpolating the received signal between two 
consecutive samples with a typical resolution between 1/4 and 1/16 of the 
sampling period. 

25 The following patents belong in this group: [Odell, "Digital beamforming and 
filtering circuit", US Patent N° 5268877, Dec, 1993], [D. Lipschutz, "Time 
multiplexed digital ultrasound beamformer", US Patent 5469851, Nov 1995], 
[Beaudin et al., "Delay generator for phased array ultrasound beamformer" US 
Patent 5522391, Jun 1996], [Wright et al. "Method and apparatus for receive 

30 beamforming system", US Patent 6042547, Mar 2000] and [M.H. Bae, 
"Ultrasonic signal focusing method and apparatus for ultrasonic imaging 
systems", US Patent 6231511, May 2001]. 

This technique, with diverse variations to physically implement the gross and 
35 fine delays, including the phase rotation by CORDIC algorithms, invariably 
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requires an interpolation stage and a selector of interpolated samples per 
channel, with the subsequent associated complexity and errors inherent to such 
systems. 

5 2. Direct acquisition of signals in the instants that correspond to the arrival of 
the signal from each focus to each element. To do this, the sampling clock 
phase is controlled for each channel with an equivalent temporal resolution 
of 1/16 and1/(64) of the fundamental period of each signal. This method 
makes the interpolation stage unnecessary, eliminating the corresponding 
10 electronic and error sources. The samples can be acquired over a FIFO 
memory from which they can be extracted to do the coherent summation 
[Chang et al., "Bandwidth sampling technique for digital focusing in array 
imaging systems", US Patent 5581036, Dec 1996]. 

15 The greatest problem associated to this technique is to achieve the sufficient 
precision in the controlled phase sampling clock generators. In [O'Donnell et al., 
"Method and apparatus for fully digital beam formation in a phased array 
coherent imaging system", US Patent N° 4809184, Feb 1989] and [Corl, "Digital 
ultrasound system with dynamic focus", US Patent 4974211, Nov 1990], m 

20 clocks derived from a master m frequency - sometimes higher than the 
sampling frequency- sampling clock are generated, each with a 2iT/m offset; a 
clock selector per channel actuated by the content of the focusing memory, 
allows associating at each moment one of the clocks to the sampling clock. The 
distribution of multiple clocks and their selection by temporal coherence is a 

25 problem, especially for the greater temporal resolutions. 

Another possibility is to store each of the channel sampling clocks in a table as 
a 1s and 0s sequence, as suggested in [J.H. Kim, T. Song, S.B. Park, 
"Pipelined Sampled-Delay focusing in ultrasound imaging systems", Ultrasonic 
30 Imaging 9, pp75-(91), 1987], but it is not a practical approach due to the amount 
of memory necessary to store all the sampling clocks of all the channels for 
multiple scanning lines. 

Another alternative calculates in real time the sampling instants by means of 
35 specialized circuits associated to each channel. This technique is described in 
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[Park et al., "Real time digital reception focusing method and apparatus 
adopting the same", US Patent 5669384, Sep 1997], [Bae, "Focusing delay 
calculation method for real-time digital focusing and apparatus adopting the 
same", US Patent 5386881, Feb 1998] and [Petrofsky, "Method and apparatus 
5 for distributed focus control with slope tracking", US Patent 5724972, Mar 
1998]. All the methods, with their variations, are based on the average point 
algorithm described in [Bae et al., "An efficient Real Time Focusing Delay 
Calculation in Ultrasonic Imaging Systems", Ultrasonic Imaging, 16, pp 231-248, 
1994]. 

10 

Although this is an interesting technique because it offers the possibility of 
integrating it in an ASIC, it has several disadvantages, such as: the need to 
program a large number of parameters for each of the scanning lines and each 
channel, the complexity of the hardware required to obtain satisfactory 
15 resolution, the lower efficiency of use of the silicon area versus that required by 
memory-based systems due to the low regularity of the hardware. It also 
requires a compromise between a high calculation speed or a reduction in the 
range of application. 

20 The electronic method and system that are the object of the present invention 
belong to this group and solve the problems described by a mixed technique 
that combines a focusing memory of a decreased size that only requires a 
fraction of bit per each focus and channel, with a high resolution of the sampling 
instant calculated in real time achieved by simple and regular circuits that are 

25 easily VLSI integrated. 
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DESCRIPTION OF THE INVENTION 
INTRODUCTION 

5 A beam former that forms the coherent composition of signals coming from a 
specific direction of propagation, received by a set of N transducers receptors or 
array is described. The methodology and apparatuses described are applicable 
to fields such as ultrasonic imaging (echography and Non-destructive Essays), 
radar, sonar, acoustics, seismology and other. 

10 

The beam former described herein operates with arrays which elements may be 
arbitrarily distributed. In particular, it achieves the dynamic focusing with 
deflection of the beam in azimuth and elevation using bidimensional arrays. 

15 The operational principle is based in sampling the signal received by each array 
element at the instants that correspond to those of the signal arrival from 
multiple foci located along a given direction of propagation. Within the limits 
established by the guiding principles of the present invention, the instant at 
which the signal that comes from each focus is received by each element in the 

20 array and acquired is determined with absolute error, in relation to the exact 
value, lower than half a period of a master clock. The master clock frequency 
can be chosen to reduce the negative incidence of the quantification lobules of 
the delays in the signal/noise ration, contrast and dynamic range of the resulting 
signal. 

25 

The sampling frequency is independent from the signal, provided the Nyquist 
criterion is verified, both for radiofrequency and for base band, as well as the 
master clock frequency. With oversampling techniques (AI conversion), the 
digitalization frequency may be the lowest frequency compatible with the 
30 desired signal-to-noise ratio. Particularly, a modified master clock can be used 
according to the method describe to do the AI sampling directly. 

Reducing the sampling frequency to the lowest possible compatible with the 
band width of the signals and/or the signal-to-noise ratio, reduces cost, 
35 energetic consumption and the volume of data to be processed. More 



particularly, with the method that is an essential part of the present invention, 
only those samples needed to do the coherent summation are acquired, without 
having to implement an interpolation process or having to acquire redundant 
information. 

5 

Also, the forming method here described, and that is also an essential part of 
the present invention, has low requirements for the focusing memory, that may 
be of a fraction of a bit per acquired sample in each channel. This high level of 
efficiency makes the method easily and conveniently applicable to various 
10 modes of operation, such as linear or sectorial scanning, with dense or 
dispersed apertures, and other applications that require multiple descriptors of 
the dynamic focusing for each configuration of the active aperture of the array. 

Said characteristics make possible the joint integration of the focusing memory, 
15 of the data acquisition, the logic for the generation of the sampling clocks and 
the integration of other auxiliary functions such as apodization, dynamic 
aperture and control of multiple channels in one sole VLSI device or 
configurable FPGA type standard, with the subsequent reduction in cost without 
losing in quality. 

20 

Description of the method 

On behalf of simplicity, the description is done for a one-dimensional array of N 
elements, not necessarily equidistantly spaced, considering an application of 
ultrasounds using the pulse-echo technique in a homogenous media with a c 
25 speed of propagation. The method described immediately below will be 
generically called Progressive Focal Correction (PFC), or the technique that 
describes the operational procedure. 

Figure 1 shows a diagram of the general geometry including the former 
30 constraints. The array of N elements numbered 1, 2, 3,..., k, ...AM and N is 
placed on the x axis with the coordinates origin arbitrarily defined, for 
convenience, in the center of the array. The purpose is to acquire a signal in the 
angular direction 9 in relation to the z axis. To do this, multiple foci are defined 
at AR intervals in the direction of the propagation (incoming signal). The method 
35 propose combines the signals coming from each foci and received by the 
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different elements of the array in such a manner that an output is obtained that 
represents the signals that would be provided by a transducer of the size of the 
array focused on each and every one of the foci. This process is called coherent 
composition of the received signals, and the result thus obtained is called 
5 dynamically focused signal with deflection. 

The AR intervals are conveniently chosen as a function of the period of a 
master clock T x =1/f x and of the speed of propagation c, such as: 

10 

(1) &R=VcT x /2 

Where v is a whole number that expresses the number of periods of the master 
clock that there are between two consecutive foci in pulse-echo. On the other 

15 hand the AR intervals must be sufficiently small to verify the Nyquist criterion: if 
T s is the temporal interval between samples and T R is the fundamental period of 
the signal, the Nyquist criterion is verified typically for T s < T R /3, considering 
the bandwidth of the habitual signals and the acquisition in radiofrequency. 
Preferably, the temporal interval between samples is a multiple of the period of 

20 the master clock, that is T s = s 7 Xj s being a whole number. 

In general, an m number of samples may be acquired between two consecutive 
foci, and choose, preferably, the value v as a function of the interval between 
samples T s = s T x for v = ms, when m and s are whole numbers. In a particular 

25 embodiment in which m =1, each sample coincides with a focus, and therefore 
all the samples are dynamically focused. When m>1, the dynamic focusing is 
strictly done on the foci, although the focus depth allows, in general, to assume 
that the intermediate samples are well focused. A variation of the method, 
called Variable Progressive Focal Correction or VPFC, varies dynamically the 

30 number of samples acquired between two foci, reinforcing the idea of 
maintaining a good dynamic focusing in all the acquired samples. 

On the other hand, and as it is well known, the value y = Tf/T x , or the 
relationship between the period of the received signal and the temporary 
35 resolution of the system, affects the dynamic range and the contrast of the 
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images constructed with multiple focused beams. In general, the value of p > 16 
(see, for instant, [D.K. Peterson, G.S. Kino, "Real-time Digital Image 
Reconstruction: A Description of Imaging Hardware and an Analysis of 
Quantization Errors", IEEE Trans. Sonics Ultrasonics, 31 pp. 337-351, Jul 
5 1984]). 

It is important that the temporal resolution as determined by T x be independent 
from the sampling period T s , that is chosen according to the Nyquist criterion, 
and therefore, it is related to the period of the signals received, T R . When the 
10 sampling is done in radiofrequency, generally Ts ~ TrI 3 or Ts= T r 14 are 
chosen to decrease the flow of data to be processed and the energy 



consumption, that is, the lowest sampling frequency compatible with the Nyquist 
criterion. In that case, when s £ 6 or s £ 4 respectively the condition of p = 16 
must be fulfilled. A particularity of the PFC and the VPFC techniques is that they 
15 allow most electronics to operate at the sampling frequency, reserving the 
higher frequency of the sampling clock for a very few elements that adjust the 
sampling instant. 

For a focus F, located on the generic polar coordinate (R, 9) and for the element 
20 k of the array located on the Cartesian coordinate (x k) 0), by the cosine theorem 



25 Where L k (R,6) represents the length of the distance traversed by the signal from 
focus F located on (R,6) to the k element. The flight time of an ultrasonic pulse 
from an emission arbitrarily chosen at the coordinate's origin to focus F and 
from the latter to the k element is: 



In case the origin of times is located on a different coordinate (for instance, 
when the space inspected is illuminated from an emitter different from that of 
the transducer receptor), the first term of the Equation (3) must be modified to 
35 represent the time of flight from the emitter to the focus F. 



(Figure 1): 




30 (3) 
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In general, the time the signal takes to reach the k element in the array from 
focus Fi is called T ki . This generic notation refers to both the geometry indicated 
in Figure 1 as to other situations in which the elements of the array are 

5 distributed arbitrarily, the propagation medium is not homogeneous or the foci 
are placed in other positions. Notwithstanding, to facilitate the description, the 
geometry of Figure 1 will be kept, the homogenous medium with a c 
propagation speed and the positions of the foci will be aligned forming a 9 angle 
in relation to the array, being placed, for the moment, at regular AR intervals 

10 from an initial distance Ro in which the firs focus Fo is placed. 

The variation on the arrival instant to element k of signals coming from two 
consecutive foci Fm and Fi, located at distances R r AR and Rj from the 
coordinates origin, respectively is dictated by the following difference: 

15 

(4) Ar tt (* ? ,0)=r t (* | ,e)-T Jt (* i -M,e) 



That, with a sufficiently small AR is: 



20 



dT (R,9) 



R = R 



(5) 
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Deriving the Equation (3) from R and substituting the value of AR obtained by 
Equation (1): 



(6) 
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This is a monotonic increasing function between R=0 and R^°° and 
independent from the speed of propagation with the following limits: 



(7) 



lim 



K 

V 1 J 



->0 



5 (8) ,??>*Kr 

Where, for a clearer notation, AT K = AT K (R,6). The value Xk/IxkI of Equation (7) 
represents the sign of x K > while G intervenes with its sign. In any case, the 
10 smallest possible value for AT K is cero, independently of the position of the k 
element. 

In this manner, the intervals between the sampling instants corresponding to the 
two consecutive foci are limited by the values expressed in Equations (7) and 
15 (8). Particularly, for 6 = 0 (there is no deflection), the function increases 
monotonically between vTJ2 and vT x when the distance of the foci varies from 0 
to * When signals are received with deflection, the range of variation is greater 
and, in a extreme case of 90° deflection is between 0 and vT x . 

20 In a possible embodiment, the positions of the foci are defined by the values 
contained in the interval [0,v], representing the incremental values given by 
Equation (6) in periods of the sampling clock. In this sense, if the whole number 
v may be expressed in b bits, each incremental delay for each foci and for each 
element may be represented with this same width, regardless of the position 

25 occupied by the array element, of the distance at which the focus is located and 
of the deflection angle. The absolute error calculated in the determination of the 
sampling instant is lower than half a period of the sampling clock when the 
result obtained by AT K (R,Q)/ T x as calculated by Equation (6) is rounded to the 
closest whole number. The disadvantaged presented by this option is that it 

30 requires a high number of bits to express the incremental delays per foci and 
per element, with the subsequent demands on memory capacity that this 
entails. 

Inversely, if there are b bits available to represent each incremental delay, the 
35 range of possible values can be chosen from two extreme situations: 



(9a) 
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(9b) 




x 



The form (9a) is useful to implement a dynamic focusing from a minimum 
distance of R min , that is, for R > R min , while the (9b) option would be preferable 
10 for dynamic focusing of the regions closest to the array, that is, for R < R max - 
There are may be intermediate ranges between these two that are more 
suitable for a particular application, where R min ^ R ^Rmax- 

However, with the Progressive Focal Correction technique, the preferred option 
15 is (9a) because it does not limit the maximum valid distance. For the regions 
closest to the array, in which is not possible to verify the inequality to the left of 
the Equation (9a), the proposed technique provides a series of alternatives that 
will be described in detail later in the paper. Therefore, if for Equation (9a): 

20 (10) a = 2*-f 

It can be described in a more compact form as: 

(11) (v-a)T x <AT K <vT x 



Equation (6) always verifies the inequality of the right side of Equation (11), 
according to Equation (8). In order to also be able to verify the left side: 



25 



30 




(12) 
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With: 
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* | Uf -xylene J 2 



Equation (12) becomes: 
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10 (14) 



1 1 



i+ 



x|cos 2 e 



V ^ 

->v - a 
2 



Where if we change the variable: 
15 (15) 



Z k ~ X k cos ^ r-R- x^senQ 
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The result obtained is: 



(16) 
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And solving for rwe obtain: 
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(17) 
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That after undoing the changes of variables done for Equation (15) leads to the 
following result: 



R>R 0 {x k ,Q) = 



v -2a 



Ivov - a 



x, COS0 
2 * 



v > 2a 



(18) 



Where l*l represents the absolute value and contemplates the most unfavorable 
case of signs relative to X k and 8. Equation (18) represents the upper limit of the 
minimum distance Ro(X k , G) from which the temporal intervals to sample the 
signal received by element k with a 6 deflection angle may be represented with 
10 a whole number of b bits with an error lower than half a period T x of a master 
clock. 

If Equation (18) is verified for all the elements of the array for all distances and 
angles of the region under inspection, the initial delay to acquire the first 
15 sample, corresponding to Focus Fo in channel k is defined by: 



(19) 
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20 

Where R 0 = Ro (X k ,G ) as defined by Equation(18). Obviously the first 
acquisition may be done at a Ri>R 0 distance, in which case the value of T 0 
(X k ,Q ) will be calculated for this new value. 

25 From this moment on, the sampling intervals between foci present a bounded 
error delimited at half a period of the master clock, that can be codified for each 
foci F/ of {Ro,0 ) coordinates in channel k as: 



30 (20) 



Ti 



Where 1*1 jj represents the rounding to the closest whole number function and 
AT Ki (Rj,G) is given by the value obtained with Equation (6). Substituting in (20) 
Equations (10) and (11): 

5 (21) Q<Q ki <2 b -\ 

That indicates the sampling intervals may be codified with b bits producing an 
error lower than half a period of the master clock. The value of Qk, is the code of 
10 focal correction that determines de number of periods of the master clock that 
the acquisition of the signals corresponding to focus / in channel k must be 
advanced. 

Acquisition of a number m > 1 of samples between foci increases the efficiency 
15 of the memory since they share one sole focal correction code of b bits between 
several samples, with a density of b/m bits per sample. If, for instance, b-1 and 
m = 4, 0'25 bits are required to focus each acquired sample. Distributing 
uniformly the advance expressed by the focal correction code Qk/ between the 
m samples is a good choice to reduce the error in sampling intervals. To do this, 
20 in a preferred embodiment the following values are calculated as follows: 




1£ j<m 



25 Where q = Q jk represents, generically, the focal correction code to simplify 
notation, and l*lti is the rounding function. The value of u(j,q) represents the 
closest whole number to the number of periods of the master clock between the 
first focus and the j sample in which j = 1, 2,..., m, achieving the most uniformly 
possible distribution of the q advance in the m samples. 

30 

The temporal interval between any given focus with the sample that follows it in 
the order j=m is u(m,q) = v -q, that is, the sample y=m is acquired in the position 
of the next focus exactly, with an advance of q master clock cycles in relation to 
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the nominal interval v between foci, according to that expressed by the focal 
correction code. 

The interval between samples, in master clock periods, is given by the following 
5 formulae: 

Au(l,q) = u(l 9 q) y = l, as per formula 
, q) = u(J , q) - u(j -l,q) 2<j<m 

(23) 

10 In the most usual case, in which the number of samples between foci is greater 
or equal to the maximum advance to be applied between foci, that is q < 2 b -1 < 
m, when uniformly distributing the q cycles between m samples, the advance gj 
assigned to the sample j can only take the values of 0 or of 1 in relation to the 
nominal interval v/m that is given for q = 0. This fact facilitates distributing the q 

15 advance uniformly between a number of samples m > q by means of a table in 
which the value g(j,q) is codified with 1 bit after having been calculated 
previously with the following formula: 

g(M)=v/m-AnC/, g ) l<j$m 9 0<q<2 b -l 

(24) 

20 To clarify the method, let's consider an example in which m=3, v =12 = 4m, b=2. 
The possible values of q are 0,1, 2, and 3 because it is codified with 2 bits, and 
v /m =4. Table I shows the results obtained when solving Equations (22) (23) 
and (24) for said values. It can be observed that advance q is uniformly 
distributed amongst the samples within the resolution of 1 period of the master 

25 clock for all the focal correction codes. The density of codification in this 
example is of 2/3 = 0.67 bits/sample. 

Table 1 

Example of uniform distribution of q advances amongst m samples 

30 
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For each value of q there are g(j]q) possible m values, each one of 1 bit, 
therefore, all the individual advances can be codified in a memory of 
5 words of 1 bit (in the previous example a 12x1 bit memory is sufficient). 

For values R < R 0 (X k , 9 ) the absolute error of the interval amongst sampling 
instants can be above the T x /2 value. In this case, the PFC and VPFC 
techniques can operate in various manners, amongst which the following 
10 strategies are most notable: 

1. Activate a dynamic aperture function in which the k element does not 
intervene until signals from a distance equal or above R 0 (X k ,Q ) are 
received. One of the possible embodiments produces for this element a 
15 number of null samples M z (k) equivalent to those produce at the output of 
the former between the beginning of acquisition T ini and the instant T 0 (X k ,6) 
given by Equation (19): 



A/ z (*) = 



20 (25) 



Where l*lt represents the function of rounding by excess to a whole number. 

It must be noted that the relationship D/ max (R 0 ), where D is the size of the 
25 aperture, determines the minimum value of F# that can be applied with the 
proposed techniques. In a practical example it is shown that F#min = 1, 0 } 6 and 
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0'4 for values of b = 1, 2 and 3 respectively, very competitive in relation to the 
values provided by other techniques. 

2. Anticipate the acquisition of the signal in the k element at a distance of R 1 < 
5 Ro (Xk,Q) chosen in such a manner that a sample in R 0 (X k) 9) is required. 

The previous sampling instants will have a certain error above the delimited 
error, but its effect on the image may be limited. To do this, the acquisition 
instant of the first sample as given by Equation (19) is then modified to: 

10 

(26) T* (^,e)^7V(^,e)^M z (i)7 51 

Where it is expressed that, for the k element Mz(K) samples are acquired prior 
to the reception with bounded error and T si is the sampling period used, which 
15 is constant in this region. The value of T s i can be equal to the output sampling 
period of the Ts beam former or, preferably, a lower value to reduce inasmuch 
as possible the errors in the sampling instant. 

3. Operate in a mixed manner, where a first range of distances [R2,Ri] is 
20 covered with dynamic aperture, another [Ri,Ro] with a limited error as in the 

previous point, and finally, a third [R 0) 0] in which the sampling errors are 
bounded within half a period of the master clock by applying the progressive 
focal correction technique. 

25 In a preferred embodiment that, in addition can be generally applied when the 
geometry of the array is different, the medium is not homogenous, (for example, 
changes in the propagation speed), or the foci are distributed at regular 
intervals but following a trajectory different from the one considered, a 
calculation algorithm based on the following steps is used: 

30 

ALGORITHM 1 

STEP 1: the times of propagation Tki of the signal corresponding to each foci Fi, 
I = 0, 1, 2,... are calculated up to each k element of the array. The 
value of Tki is obtained in each case when considering the geometry of 
35 the system and the speed(s) of propagation. 



20 
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STEP 2: For each k element the U k i intervals between the instants of arrival 
of the signals are calculated for each of the two consecutive foci 
expressed in periods of the master clock, such as: 



(27) 



T -T 
ki k,i-\ 



'X 



i>1 



STEP 3: The value v, or nominal interval, in periods of the master clock 
between foci is chosen from the results of the previous equation as: 



v>maxfj/ w ) VJt.i 



(28) 



STEP 4: If the number of bits available to codify the focal corrections is b, 
determine for each k channel the index i=h k from which: 



15 (29) 



v-£/ w ^2°-l 



i*h, 



STEP 5: For foci F, with /< h ki> the number of samples Mz(k) to operate with 
dynamic aperture or advance acquisition in channel k is: 



20 (30) 



STEP 6: The focal correction codes for each k element corresponding to focus 
i> hk are calculated as: 



25 



(31) 



G*=v-t/« \>h k 



30 



Expressing the number of master clock cycles with which the acquisition must 
be advanced for the signal corresponding to focus / in channel k. 

In another preferred embodiment, during the acquisition of the signal, the 
distance between two consecutive foci Fm and F, is varied at moments selected 
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by increasing the number of samples acquired between them. This variation is 
called Variable Progressive Focal Correction technique or VPFC in which the 
number of samples acquired between two consecutive foci goes from being a 
constant m to a be a function of the number of current focus, that is M = m(i). 
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Equation 18 indicates that the value of R 0 grows along with the value of v. On 
the other hand, the minimum value of K is determined by Equation 1, when one 
sample is acquired for focus, with a temporal interval between samples of Ts = 
sTx. Initiating then the acquisition with M = m (1) = 1 from a distance R 0 (1) as 
determined by the value v = s to which for all the elements: 



j-2a 



=x, cosG 
2 * 



I k j 



(32) 



15 Where a = 2* - 7. The acquisition continues maintaining M = m(1)=1 until 
certain focus is reached at a distance R= Ro (2) based on which the following is 
verified for all the elements: 



20 



* 0 (2)£ 



(33) 



25 -2a 



\4los-a 2 



jc^ cosG 



V* 



starting on a section in which 2 samples may be acquired per focus, that is M 
=m (2) = 2, while maintaining the bounded sample error. After a certain time 
another focus, R= R 0 (3) will be reached and can be operated with M=m(3)=3, 
25 and so on, in such a manner that, generically, from a focus located in R= R 0 (m) 
from which m samples can be acquired verifying that: 



(34) 



ms-2a 



l4ams-i 



= Xy cosG 

2 * 



V* 
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In summary, with the Variable Progressive Focal Correction Technique the 
number of samples acquired between foci is increased with time, and the focal 
correction codes can continue to be expressed with b bits. It must be noted that 
5 the condition (34) must be verified by all the elements for each value of m. 

In a preferred embodiment a 1 bit field called focal offset code J h is associated 
to each focus F,. When J,=0 the current value m(i)=m(i-1) is maintained, and 
when Jj=1 it is incremented as: m(i)=m(i-1)+1. Algorithm 2 provides a possible 
10 method to calculate the values for J,. 



ALGORITHM 2 

STEP 1: based on the minimum distance established by the application, Rmin, 
the larger value m-m 0 is determined and serves to verify Equation (34). 

15 

STEP 2: Initialization: m=m 0 . Auxiliary variables: f=1 (current focus number), i=1 
(current number of samples). 

STEP 3: Repeat: 

20 a) Calculate R 0 (m) given by Equation (34) and the value of AR = cTs/2 

when c = speed of propagation, Ts = sampling period. 



b) If Rmin + (i + m)AR> (m+1) R 0 

{Jf= 1; m=m+1 } 
25 In another case {Jf = 0} 

c) f=f+f; 



d) /=/+m; 

until i>A/s, where Ns = number of samples to be acquired. 

The number of focal offset codes is equal to that of foci for any acquisition, 
which can be easily done by associating to each of the focal correction codes 
Qki the appropriate focal offset code J k i = Ji in each k channel. 
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This scheme is appropriate since, in general, the signals corresponding to each 
focus arrive at different moments to the elements and it is easy to process. In 
this manner, the focal correction memory has a word band of b+1 bits 
distributed in two fields; b bits for the Q ki code and 1 bit for the J, code. 

5 

An alternative that provides for greater efficiency, stores in a separate memory 
the focal shift codes J,, which content is shared by a set of channels. To solve 
the problem of non-simultaneity in the foci acquisition each channel has a small 
tail where the J { codes obtained automatically from the memory where they are 
10 stored are entered. Each channel uses the value of J, at the head of the queue 
each time it reaches a new focus. In this manner, one sole central memory 
stores one J, code per focus and some distributed small FIFOs facilitate the 
assignation of these common codes to each channel based on need, as it will 
be described in greater detail later in the system description section. 

15 

On the other hand, and since the distance between foci may be modified 
dynamically, for most of the practical applications one encoding of the focal 
correction with b = 1 bit is sufficient, and this provides a high efficiency in the 
use of the memory and a simplification of the correction circuits. In particular, it 
20 is necessary to immediately do the uniform distribution of 0 or 1 period of the 
master clock sample between m samples. 

In view of the description given above, the Progressive Focal Correction 
Technique provides, with its Variable option, a method that effects the dynamic 

25 beam forming at reception with a high temporal resolution within a great range 
of distances. Additionally, it lacks acquisition redundancies, it does not require 
processes for sample interpolation, it allows for maintaining the lowest sampling 
rate compatible with the Nyquist criterion, and, therefore has a low energetic 
consumption. Its realization is compact, due to both an adequate balance 

30 between the calculation in real time of the corrections and the content of 
memories pre-calculated in non real time, and to the very efficient encoding of 
this information. 
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All these characteristics lend themselves to configure a versatile beam forming 
system at reception with dynamic focusing, deflection, dynamic aperture and 
apodization which constructive details are detailed below. 

5 Description of the system 

The system to which the present invention refers is based on the Progressive 
Focal Correction technique described above, that may or may not include the 
Variable option, and that does the dynamic focusing deflecting the signals 
received through a set of transducers or array, including the apodization and 
10 dynamic aperture functions. The system will be described in sufficient detail to 
be understood and reproduced by personnel with knowledge of digital 
electronics. 

In reference to Figure 2, a system of beamforming in reception is composed, in 
15 the most general case, of an array (10) composed of N elemental transducers 
(11) or, simply elements that receive outside signals numbered 1, 2,...k,... f N. 
Frequently, these elements are, also, emitters of signals to form a beam at 
emission, and passing to be receivers once this function is completed. The 
elements necessary to form the beam at emission have been omitted in the 
20 figure. 

The signals received by the various elements are conditioned and amplified by 
the AAS devices (12), that deliver the analogic N signals a 1t a 2 , ... a k} ...,a N to the 
corresponding analogic-digital converters A/D (13). The AAS devices usually 
25 have, in addition to the amplifying function, other function that include 
compensation of attenuation with the distance by means of a gaining-time 
control, filtering, and on occasion the unfolding of the components in phase and 
cuadrature. In this last case each analogic signal ak is formed by the pair (i k , 
that represent, respectively the components in phase and in cuadrature. 

30 

Once conditioned, the analogic signals are digitalized by the A/D converters. In 
the case of signals in base band with pairs (i k) qd two A/D converters may be 
used per signal or they can be multiplexed in one sole A/D converter between 
two signals, the even samples for the i k component and the odd samples for the 
35 q k or vice versa. 
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Another alternative is using AI converters, and represent the analogic signals 
a k with a continuous bit sequence. In this particular case the Progressive Focal 
Correction technique and its variations are especially useful for adjusting the 
5 sampling instants to the arrival signal instants to each element with the 
resolution of a master clock which period is adapted, cycle by cycle, to that 
required by such AI converters to obtain a given signal to noise ratio. 

The techniques mentioned are known by the experts in the field and are also of 
10 general application. The differences reside mainly in the CONF device (14) or 
beamformer, which function is to combine the N numerical sequences e1, 

e2,....ek e N provided by the A/D converters to create one sole exit sequence 

r that represent, in digital form the signal received by the array and focused, 
ideally, in all the points of one given direction. 

15 

Each sample delivered by the A/D converters is a word of w bits, typically 6 < w 
<12 for the instant sampling converters and 1 £ w <2 for the AX converters, 
although the system can operate with other word widths. 

20 It is precisely the manner in which the CONF device performs the function of 
combining coherently the N input entries ek to generate the output r what 
differentiates the present invention. 

The operation of the CONF device (14), which is the subject of this patent, is 
25 based on the Progressive Focal Correction technique, which is one of its key 
aspects that contributes several advantages versus other methods. 

Therefore, according to other methods described in the Prior Art section, each 
sample of the numeric sequence e k provided by the /c-teenth A/D converter is 
30 subjected to a delay process as shown in Figure 3. Each element is associated 
to an individual processing channel (21). The content of a dynamic focusing 
memory (25) is used by a control device (24) to generate signals that produce a 
gross (22) and a fine (23) delay. The gross delay has, typically, the resolution of 
one period of the sampling clock, and the fine delay is a fraction of that. These 
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delays are obtained with various methods that have been described in the 
specialized literature and are well known. 

The content of the dynamic focusing memory (25) is such that, for each sample 
5 and for each channel the delays introduced align in time the outputs f^.-.f/v 
that correspond to the successive foci F 1} F 2 ,....F N . These signals are summed 
in (20) to obtain the dynamically focused r sequence. 

The greatest advantage provided by these techniques is that only on sampling 
10 clock, common to all the channels, is needed. However, it suffers from various 
disadvantages derived from the need to having to carry out interpolation 
processes in order to obtain sufficient temporal resolution in the delays 
introduced in the signals, with the subsequent increase of electronic complexity 
and the possibility of introducing errors. Also, they require a considerable 
15 amount of memory to represent the gross and the fine delays. 

The system that is related in this present invention operates in a different 
manner, by sampling the signals directly in the instant they arrive at every one 
of the array elements, coming from each of the foci located along a given 
20 direction. Because the introduction of delays to the signals is avoided the 
interpolation processes are unnecessary, which serves to avoid possible errors 
and eliminate the electronic circuits necessary for that function. Also, only those 
samples that are necessary to obtain a coherent composition are acquired, 
eliminating redundancies. 

25 

The initial cost of the method subject of this patent resides in that each A/D 
converter must operate with a different sampling clock of a non-uniform 
frequency. As shown in the description below, obtaining this characteristic is a 
simple matter when using the methodology described as Progressive Focal 
30 Correction with the Variable option, and it can be integrated in multiple channels 
in one sole device of the FPGA standard type and commercially available. 

Referring again to Figure 2, the proposed method requires that the 
beamforming device CONF (14) provides an independent and of a non-uniform 
35 frequency sample clock CK 1} CK 2 ...DK k , CK N for each A/D converter (13) that 
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enables the acquisition of analog signals a*, a*... a*,.... aw provided by the AAS 
signal conditioner (12) in the instants that correspond to those of the arrival from 
each foci to each element (11) of the array (10). The samples obtained by the N 
A/D converters (13) are processed to obtain their coherent sum, that is, an 
5 output value r that is the sum of the N samples corresponding to the signal 
received from one same focus by the N elements of the array. 

In the present invention, this methodology is applied in modular form by 
associating a processing channel to each of the array's elements. Figure 4 
10 shows the structure of one of the channels (30) that processes the signal 
received by the k element (11) of the array. All the channels have a structure 
identical to that shown. 

As in other cases, the signal is amplified and conditioned by the AAS device 
15 (12) by known methods that include electronic devices to amplify, filter and 
manipulate analogically the signal. The signal a k , digitalized by the A/D 
converter (13), is obtained at the AAS exit. The samples are acquired at the 
instants defined by a clock generator (70) that produces a sequence of pulses 
CK k not uniformly spaced, which detailed description will be provided later. The 
20 clock generator is enabled by the HFD signal generated in the HFOC device 
(80) that marks the instants at with the channel begins and ends the acquisition 
of the dynamically focused signal. 

The function of the clock generator is to guarantee that the sampling instants 
25 determined by CKk correspond to those of the arrival of the signal from each of 
the foci to the k element, and operates according to the principles described 
above with the name of Progressive Focal Correction Technique and its 
Variable option. 

30 With this, the sequence of e k samples obtained at the exit of the A/D converter 
corresponds precisely with the ordered values of the signal received by element 
k coming from each of the focus. These values are processed by an apodization 
block and dynamic aperture APD (60) that allow the enhancing of the beam 
characteristics and that produce the f k sequence. Because there is a 1:1 

35 correspondence in the e k to the f k samples, this sequence is also ordered with a 
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reference to the arrival instants of the signals to element k from each focus. The 
successive values of f k are transitorily stored in a tail or FIFO memory that 
absorbs the difference in the instants in which the samples corresponding to a 
same focus in different channels. 

5 

The key element is then the sampling clock generator (70) that produces the 
digitalization signal CK k of the k-teenth A/D converter so the samples obtained 
correspond to the values of the signal received by the k element from each of 
the foci. Figure 5 shows the beginning diagram of one of the preferred 
10 embodiments of the clock generator CK ki which function is based on the 
method previously described, the Progressive Focal Correction technique for 
those cases in which the focal corrections are codified with b-1 bit. 

The circuit shown in Figure 5 only contains the basic elements that are 
15 necessary for its correct operation, the auxiliary devices having been omitted to 
make the description more clear. Said circuit has a REG-A registry (71a) in 
which a value related to the v parameter or nominal interval between foci is 
loaded in periods of the master clock. The output V has multiple lines that 
determine the number of stages that make up the shift register of programmable 
20 length SHR-A (72). This device has an exit p that is a delayed version of a 
quantity of V periods of the c clock of its input CK k . This and other devices in the 
scheme shown operate synchronically with the master clock c of period T x . 

The input signal HFD enables the operation of shift register (72) and of the flip- 
25 flop (77), in such a manner that these devices keep their actual status 
regardless of the value of their input when HFD = 0. When HFD =1 a dynamic 
focusing is enabled and then interrupted when HFD returns to zero. This signal 
is generated in the HFOC device (80), Figure 4. 

30 The signal P is also delayed 1 period of master clock by the flip-flop FFA (73) 
producing the pr signal. The MUX (74) multiplexer provides an output s that is 
either the signal p produced by the shift register or its delayed version pr, based 
on the status of the selection line Q: 

35 
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(35) s=p if Q=0; s=pr if Q= 1 

The status of the Q selector is determined by the contents of a MEM memory 
(75) in the F direction provided by the CNT counter (76). In the present case, k 
5 is a simple line (1 bit). 

The memory is previously loaded, by circuits that are not shown, with the focal 
correction codes Q ki , where k is the number of the current channel and / is the 
ordinal number of the focus, with i>1. The circuits for the signal acquisition at 
10 the first focus (i = 0) that will be described later, determine the activation of the 
HFD signal. 

The output of the multiplexer (74) is registered in the flip-flop FFB (77), which 
output C K k is the sampling clock for the A/D converter in the current channel k 
15 (13), shown in Figure 4, that actuates the CNT counter (76) and is the input for 
the shift register SHR-A (72). 

In order to operate, the shift register SHA-A (72) is initialized at '0', the flip-flop 
FFA (73) at '0\ the flip-flop FFB (77) at 'rand the counter (76) is initialized with 
20 the direction of the memory where the focal correction code corresponding to 
focus Fi (i=1) is stored. The auxiliary circuits that needed to establish these 
initial statuses are omitted in the figure. Therefore, the initial values of the 
signals are: CK^~1 } p=0, pr=0 and s=0 (independently of the status of Q). 

25 From this initial state, CK k =0 is loaded on the cycle of the master clock c that 
follows activation of the HFD-1 signal, since the input of the flip-flop FFB (77) is 
in that status and the clock enabled for its input ce. During the first V periods of 
the clock c the signals p=0 and pr=0 are maintained. In the following cycle, p=1 
is entered, reflecting that the initial value of CK k =1 has propagated by the shift 

30 register (72), but is still maintained at pr=0. In this same cycle, if Q=0, s=pr=0, 
but if Q=1, s=P=1, therefore, in the following cycle, the exist CK k =1 only if Q=1; 
pr=1 must also be included in this cycle and therefore, if Q=0 the CK^-1 output 
must be entered in the next cycle. Once that CK k =f the process repeats itself. 
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That is why the number of n cycles that pass from the previous value CK k =1 
and the next value are determined by the value o Q in the following manner: 

(36) lfQ=0: n=V+1 

5 lfQ=1: n=V+2 

According to Equation (20), when Q=0 the number of cycles between foci is 
n-v, and when Q=1 it must be in view of which: 

10 (37) V=v-2 

Which indicates the initial value to be loaded in the REG-A registry (71a) as a 
function of the nominal distance between foci v and expressed in periods of the 
master clock. 

15 

Each time that a sampling clock (CK k =1) is produced, the F direction provided 
by the CNT (76) counter is increased and the memory MEM (75) places at the 
output the code corresponding to the focal correction to be applied to the next 
focus. 

20 

It is observed that the set FFA (73)-MUX 874)-FFB (77) configures a shift 
register of variable length according to the code Q, that may be used as the 
base for the actions describe below, although there are more efficient 
alternatives. Figure 6 shows a scheme of another preferred embodiment that 
25 does not limit the number of bits with the focal correction codes Q represented. 
In this case the previous FFA -MUX -FFB set has been substituted by a shift 
register SHRB (78) which programmable length is determined by the 
complement to the focal correction code Q as obtained by the inverters (75a) 
that produce the outputs A/Q. 

30 

In this case the interval between the sampling instants in two consecutive foci in 
the current channel are determined by the multibit value of Q as delivered by 
the memory (75). Here the shift register SHRB (78) is initialized as T. 
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The operation is, otherwise, similar to that of the previous case except that the 
number of cycles is determined by the V and NQ values according to: 

(38) n=V+NQ 

5 

Taking into account that A/Q=2* - 1 -Q, the initial value to load in the REG-A 
register (71a) as a function of the nominal interval between foci in periods of the 
master clock is: 

10 (39) V=v-2 b +1 

In a practical embodiment, the inverters (75a) are unnecessary if the 
complementary values of NQ instead of the values of Q are stored in the 
memory of focal correction codes (75), maintained in the description to lend 
15 coherence to the notation. Also, in the practical embodiments of the circuits 
showed in Figures 5 and 6 the auxiliary logic that avoids situations in which all 
the outputs are at '0' -and would prevent the output of circuits- is also included. 

The logic described allows the focusing of each of the acquired samples. When 
20 a periodical focal correction is preferred, sharing the same focal correction code 
Q each m samples, the previous circuits are modified in the manner shown in 
Figure 7. The value of m is loaded in an RCM register (71a) and the value 
stored in M is used by a control logic CTRL (90) together with the Q value of the 
focal correction code delivered by the memory (75) to generate the ng and the 
25 cef signals. The first controls the length of the shift register SHRB (78) making it 
equal to ng+f, while cef enables the CNT (76) counter, which provided the F 
direction to the focal correction memory (75). 

The shift register SHRB (72) has a programmable length based on the V value 
30 programmed in the REG-A (71a) register, with a number of stages equal to 
V+1. The length of the shift register SHRB (78) is also programmable with a 
number of stages equal to ng+1. 
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The number of master clock periods that occur between two consecutive CK k 
cycles is equal to the number of interposed stages between the input of SHRA 
(72) and the output of SHRB (78) that according to that explained above is: 

5 (40) n=V + ng+2 

Programming in the REG. A (71a) register a value of: 

(41) V=v/m-3 

10 

Where v/m represents the nominal interval between samples in master clock 
periods, and the resulting Equation (40) is then as follows: 

(42) n =v/m + ng-1 

The logic of the CTRL (90) control shown in detail in Figure 7, controls the n 
interval by means of the ng signal, distributing uniformly the advance of the 
periods of the master clock Q between the m samples that are acquired 
between two consecutive foci. 

To do that, the CM content of a module m counter CMOD (91) that is updated 
with each CK k clock by its enabling input ce, provides part of the DM direction of 
the MDM memory (92). The other part of the direction is the focal correction 
code Q. The direction provided to the MDM memory is: 

(43) DM(j,Q) = (j-1) 2* +Q 

Where }=CM+1 indicates the order number of the current sample. Programming 
in non real time the contents of the DM(j,Q) directions of the MDM memory with 
30 the g(j,Q) value that are obtained after applying Equations (22) to (24), the g 
output represents the individual advance to be applied to each of the samples 
during acquisition. The inverter (93) obtains the ng signal, complementary of the 
g signal so according to Equation (42): 

35 IF g=0 — > ng-1 — ► n=v/m 



15 



20 



25 
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(44) IF g=1 -> ng=0 -> n-v/m-1 

producing a 1 clock cycle advance when g(j,Q)=1 and the nominal interval v/m 
is maintained when g(j,Q)=0, as described in the Progressive Focal Correction 
5 method. 

Each time that the CMOD (91) counter overloads (that is, it switches the value 
m-1 to 0), a cef pulse is produced that enables the increment of the foci counter 
CNT (76) providing a new F+1 direction to the MEM memory (75). In this 
10 manner, the process is repeated with a new Q value that has been provided by 
the MEM memory (75) with the content of a new direction. 

However, additional observations must be made. The first observation is that, 
evidently, the inverters (93) are eliminated if the values complementary to the 
15 g(j,Q) values as obtained by the Equation (24) are propagated. Another 
observation is that the MDM (92) memory has a small size; for instance a 256x1 
bits memory may be appropriate for values m< (64) with b=2 or for m< 32 with 
b=3. 

20 The variant called Variable Progressive Focal Correction is derived from this 
scheme. To obtain it, the RCM register (71b) is substituted by a counter so that 
the number of m samples amongst foci can be increased in time starting from a 
minimum initial value ml that can be, typically, equal to 1. 

25 As described before, and according to Equation (18), the minimum distance R 0 
from which the absolute error of the sampling instant is kept bounded within +/- 
T x /2 and grows with the value of the interval between foci v, expressed in T x 
periods of the master clock. Reducing the distance between foci in the region 
closest to the array my programming a low value for mi decreases the value of 

30 R 0 . The instants at which the RCM counter (71b) is incremented are determined 
by the instants at which the signals coming from foci located at a sufficient 
distance that one more sample can be acquired between the foci at all the 
channels. The ordinal number of the foci in which this occurs is calculated in 
non real time and is coded in a table that is common to all the channels, to a 
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subset of channels or to an individual channel, establishing a balance between 
efficiency of the use of memory and modularity. 

In a preferred embodiment, and in acquisition time, the foci in which the RCM 
5 (71b) counter has been incremented, are determined by the value of a focal 
offset code J\ expressed with 1 bit, where / = 1, 2,..., associated to each foci: 
the actual value of m contained in the RCM counter (71b) is incremented if J { =1, 
it maintains its current value if Ji = 0 until the next focus, where the process is 
repeated. The RCM increment is only enabled, then, when the sampling clock is 
10 produced to acquire the signal corresponding to a focus. 

In a particular embodiment, the focal correction memory (75) provides two 
fields: one J k i =J\ of one bit indicates whether the RCM counter (71b) is to be 
incremented at this focus and the other, Q ki of b bits, indicates the advance to 

15 be distribute amongst the current m samples. The figure shows the J = J\ signal 
with a broken line. In this case, preferably, focal correction codes Qki with b-1 
bit is used, since the distance between foci can be controlled by controlling the 
number of interleaved m samples. In that case, the distribution of the advances 
expressed by Q k i are of 0 o 1 master clock cycles, and are done in the sample 

20 which order is half of the content of the RCM (71b), although other simpler 
alternatives that have little influence on the image's quality are also possible. 

Nevertheless, and to improve the efficiency, another alternative embodiment 
also shown in Figure 7, the focal offset codes J\ come from a separate memory 

25 MJ (59) -see Figure (10), which output is shared by multiple channels by the 
input J to a local FIFO queue (79). During the acquisition time, a reading of a 
new focal shift code Ji is read into the global MJ (59) memory shown in Figure 
(10), when some of the channels that services detects that the local FIFO (79) 
is empty. The supplied Ji code is simultaneously written in all the local FIFOs 

30 (79) to be used by each channel when they acquire the signal of a new focus, 
that is, their cef signal is activated. It is at this time that the value at the head of 
the line Jc is used to increment or not the m value contained in the RCM 
counter (71b). 
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It must be noted that the local FIFO queues (79) are small, typically less than 16 
x 1 bits, because they only need to store a number of 1 bit codes equivalent to 
the maximum difference of foci acquired in different channels within a same 
module. This alternative consumes fewer resources than the previous one, in 
5 which the focal offset codes Ji were stored in the focal correction memory (75) 
of each of the channels. 

With the Progressive and the Variable Progressive Focal Correction Technique 
the use of focal correction codes is made more efficient, distributing it between 
10 a small number of samples at the beginning of the acquisition (at the limit, 1 bit 
per sample) and between a greater number of samples as the distance of the 
foci to the array increases (a decreasing fraction of 1 bit per sample). 

On the other hand, it is necessary to establish with precision the instant at 
15 which the samples that would correspond to the first focus - with coordinates 
(Ro,6) would be acquired, and based on which the focal correction operates 
according to the described methodology. This instant is determined by Equation 
(19) for the geometry shown in Figure 1. It must be mentioned that the initial 
focus may be located at a greater distance based on the needs of the 
20 application. Also, the use or not of a variable number of samples between foci 
will be considered to establish this initial value. 

This function is done in the HFOC (80) ) of Figure 4, which will be shown in 
greater detail in Figure 8. One embodiment uses, simply, a counter actuated by 
25 the master clock that is initially loaded with the value N A (k). 



Where T 0 (k) represents the interval between the origin of times until the signal 
that corresponds to the F 0 focus arrives to the k element, as determined by the 
30 geometry and the speeds of propagation, or as expressed by Equation (19) for 
the configuration shown in Figure 1. 




(45) 
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The signals that correspond to the first F 0 focus are acquired in channel k when 
this counter is overloaded, and enabling from this moment on the actuation of 
the progressive focal correction circuits. The absolute error of the acquisition 
instant at this first focus is lower than half a period of the master clock in all the 
5 channels, however, to maintain this temporal resolution this counter must be 
actuated at the master clock's frequency, which is normally high and 
considerably above the sampling frequency. 

To avoid this inconvenience, a preferred embodiment decomposes the counter 
10 of N A cycles in two parts, as shown in Figure 8, in such a manner that (channel 
index k omitted) is as follows: 

(46) Na - k Nai + Na2 + 1 

15 where K, N A i and N A 2 are whole numbers. The value of N A i is loaded on the 
CNT-A (82) counter by means of the ida signal. When it is enabled by actuating 
the cea signal, the counter is actuated with a cb clock which frequency is K 
times less than that of the master clock, c, and thus the mentioned difficulties 
are avoided. 

20 

After N A i cb cycles equivalent to KN A i periods T Xi the counter reaches the 
terminal value that lifts the signal tea = 1. In the next cycle of the cb clock the 
type RS (83) flip-flop FF is activated, and the signal aj=1 is placed. This value is 
propagated throughout the shift register SHRG (84) which length establishes 

25 the value programmed for N A 2, as actuated by the master clock c. After N A 2 + 1 
periods of the master clock have passed, the aj = 1 value reaches the output of 
the shift register, and the output goes to HFD = 1, enabling the dynamic 
focusing circuits after N A periods of the master clock have passed, and as 
expressed by Equation (46). In this manner, only the shift register SHRG (84) 

30 must operate at the high frequency of the master clock. 

This state remains stable until an N F number of focal corrections have been 
done as programmed in the CNT-F counter (86). Each focal correction is 
indicated by an input signal cef } that provides the control logic CTRL (90) shown 
35 in Figure 7 when multiple samples are taken between foci, or it is the sampling 
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clock CK k if the signal is only acquired between foci (Figures 5 and 6). When 
CNT-F reaches the terminal count, the tcf-1 signals is raised, provoking the 
setting to zero of the flip-flop FF 84 through its R input. When the function AND 
(81) is completed, the HFD signal returns to 0, ending the process of 
5 progressive focal correction. 

However, the acquisition of signals may continue with a constant sampling 
period also equal to the nominal v/m value until the number of samples N s to be 
acquired has been obtained as programmed in the CNT-S (88) counter. The 
10 overflow of this value activates the FIN signal that terminates completely the 
acquisition of the signal in this channel. 

It is important to observe that once the progressive focal correction process has 
been completed when HFD=0, the F direction of the focal correction codes 
15 memory (75) becomes frozen at the last value reached by the counter (76) as 
shown in Figures 5, 6 and 7, and thus making possible that the next acquisition 
set to the focal correction values stored next can be effected without the need to 
program a new initial value for F. In this manner it is possible to operate with 
various strategies, for example: 

20 

a) Generation of sectorial image where each acquisition is done with a 
different angle of deflection. The focal correction codes that correspond to 
the successive deflection angles are stored sequentially in the memory, 
avoiding the need to program the initial direction F of the focal correction 

25 codes memory before initiating each acquisition. Contrary to other 

techiques, this technique allows maintaining the bounded and the constant 
sampling error for all the deflection angles and at all the distances. 

b) Linear scan and other manners of acquisition where an active sub- 
30 aperture composed by a subset of the elements of the array changes in 

each acquisition eliminating certain elements and introducing other. 
Similarly, the focal correction codes that correspond to the successive 
configurations are sequentially stored in the memory which avoids having 
to load the initial direction F of the focal correction codes memory before 
35 initiating each acquisition. 
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c) Mixed strategies in which the signal coming from multiple directions is 
acquired using the same focal correction codes and different values N A (k) 
for different directions to those used to acquire the first focus. In this case, 
5 the value of the initial direction must be re-established before the 

acquisition. 

The APD device (60) shown in Figure 4 does the apodization and dynamic 
aperture functions. Apodization, which consist in multiplying the signal received 

10 in each channel by a positive real number A k from 0 to 1 improves the 
characteristics of the formed beam and the contrast of the image because it 
reduces the lateral lobules. The dynamic aperture makes possible maintaining 
the constant lateral resolution operating with a constant F# number that requires 
the size of the active aperture increases progressively in time. The dynamic 

15 aperture also facilitates the application of the Progressive Focal Correction 
technique in the regions closer to the transducer. 

Figure 9 shows a block diagram of the APD device (60) that uses a multiplier 
MULT (61) for the apodization function with which it obtains the product of each 

20 of the values of the e k sequence times the A k coefficient of apodization assigned 
to the k channel. In a preferred embodiment this coefficient is a whole number 
without sign of a bits. If the e k values are expressed with w bits, the most 
significant w bits of the a+ w bits that make up the product are maintained and 
rounded as a function of the value of the first ignored bit and the sign of e k} 

25 including the corresponding logic in the MULT (61) device itself to obtain the h k 
output. The function is then: 

(47) 




30 It must be noted that the apodization function may be enabled or not with simple 
control circuits that are not shown in Figure 9. Basically, output h k is input e k if 
the apodization function is enabled, in other case it responds to Equation (47). 
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This function is done with a simple multiplexer governed by the apodization 
enabling signal. 

For the dynamic aperture function there is a binary counter CNZ (62) ) that is 
5 initially loaded with value MZ(k) or number of null samples as expressed by 
Equation (25) for channel k. This counter is actuated with a clock with a nominal 
sampling frequency, cs and is enabled by activating the HFD-1 signal that 
initiates the dynamic focusing process and is provided by the HFOC (80) device 
shown in Figures 4 and 8. The signal HFD=1 acts on the S signal of the FF flip- 
10 flop (64) sending the signal z=1 that selects the null input (63) to the multiplexer 
MUX (65). 

The acquisition and processing of signals in all channels is initiated in this state. 
In those that the CNZ counter (63) has not reached the terminal count {tc 

15 signal), the output samples in fk will have a null value, and therefore the system 
operates as if the associated element was not active. After MZ cycles of the 
sampling clock, during which other null samples over fk have been produced, 
the fc=1 signal is raise setting to zero the FF flip-flop (64), and enabling the 
output of the h k values obtained by the multiplier (61). In this manner a 

20 programmable number of null samples are provided at the beginning of each 
acquisition in each channel, which permits programming the characteristics of 
the dynamic aperture. The ft sequence is then: 

f k = {0, 0, ...,0, AtficfNZ+l), A& k (NZ+2) 9 ...} 

(45) 

25 Where e k (i) is the sample i of the e k sequence and therefore it corresponds with 
the / sample of any other secondary sequence. In this manner the order is 
maintained to later do the coherent composition of all the signals, null or 
otherwise. 

30 Another important aspect of the present invention is its capacity for modular 
embodiment, a concept that is shown in Figure (10). A preferred embodiment 
groups in a module (50) all the processing functions of the signal for a subset of 
PS elements of the array, as well as other functions of global control. In 
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particular, a module that includes the MJ memory (59) that provides the focal 
shift codes J together with the circuits necessary for the reading and initial 
loading. Also, each module (50) contains S sub-modules (40) where each sub- 
module process the signals from P elements and a combinator (55) that allows 
5 chainlinking multiple modules. 

Internally, the S sub-modules (40) that make up one module are chained to 
obtain a coherent composition of the PS channels in cascade through the 
FIFOs (52). This structure facilitates the operation in parallel of all the channels 
10 of digital processing where the FIFOs compensate automatically the differences 
in the instant of acquisition of the signals in the different channels. 

To do this, each of the elements (11) of the array (10) provides the signal 
received to one of the P elementary processors PE (30) integrated in a sub- 
15 module (40). Each elementary processor PE effects conditions and digitalizes 
the signal at the arrival instant of the signal to its corresponding element, 
following the methodology described above, and includes the devices shown in 
detail in Figures 4 to 8. 

20 Each PE provides at its output an ordered sequence of values that correspond 
to the signals received by the element associated from each of the foci or 
interleaved samples. These outputs are added (41) to obtain the coherent 
composition of the P signals that correspond to one sub-module. 

25 Device (42) adds up the results obtained locally at the sub-module with those 
already placed in its FIFO queue (52) by its predecessor. In this manner an 
ordered sequence that corresponds to the focused signals of its P elements 
plus those of the sub-modules that precede it in the chain is obtained at the 
output of each sub-module in the order of signal arrival from each foci. The first 

30 sub-module in the chain receives a value of '0' (51) in substitution of the FIFO 
(52), since there is no preceding sub-module. 

Also, the combiner (55) add in (58) the results obtained by this module and that 
are stored in the FIFO (56) with which the preceding module has placed at the 
35 FIFO's queue (57). Since in both cases the results are ordered sequences of 
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the values corresponding to the successive foci, the output sequence of the 
adder reflects this same order. The chainlinking of L modules (50), each with S 
sub-modules (40) with P elementary processors in each sub-module provides a 
system of dynamic focusing with the Focal Correction technique Distributed for 
5 N=LSP elements. 

It is therefore a structure that presents several advantages: 

■ It is modular, which facilitates the development of circuits and endows it 
1 0 with greater flexibility of application, 

■ It process in parallel, which allows it to attain high rates of calculation, 

■ It has an unlimited capacity for expansion, logically-wise. 

Figure (11) shows one of the possible embodiments of the adder (41) integrated 
15 in one of the sub-modules (40). This adder receives the P sequences ei a e p 
corresponding to the values of the samples obtained by the P A/D converters 
corresponds to the first focus, the following to the second focus and son on. 
This reasoning is extensive when applied to samples obtained between foci, 
with or without the option of the Variable Progressive Focal Correction method, 
20 however, the sampling instants are not simultaneous in all the channels but it is 
possible that one channel acquires several successive samples while other 
channels wait for the appropriate signals to arrive at the associated element. 

This phase offset is automatically corrected with the distributed FIFO memories 
25 (43), (45) and (47) together with the adders (44), (46) and (48). FIFO (49) 
performs the same function but in reference to the sample sequences of two 
sub-modules. 

The FIFOs are associated in pairs which outputs are added when both have, at 
30 the least, one sample stored. This situation is detected with a simple control 
circuit that executes a simultaneous reading in both FIFOs when they are not 
empty and writes the result obtained by the corresponding adder in the FIFO 
connected to its output Each reading liberates one position of memory in the 
FIFOs that may be occupied with the writing of a new sample. Each writing 
35 action occupies a new position in the memory of the appropriate FIFO. 



42 



Readings in the associated FIFOs are simultaneous but not so the writings that 
depend on the sampling instants at the respective A/D converters, that are also 
derived from the arrival instant of the signal from each foci to the element. 

5 The capacity H of each FIFO (43) must be sufficient to store the samples 
acquired until there is a sample in the associated channel, at which time the 
adder (44) begins operating, releasing the position occupied by the 
corresponding sample and transferring the result of the addition operation to the 
FIFO (45). 

10 

Similarly, the FIFOs (45) must have a capacity of 2H, FIFOs (47) a capacity of 
4H and so on. In this manner the amount of memory distributed as FIFOs that is 
required by a sub-module with P inputs is: 

15 

H>P+2HP/2+4H*P/4+..=HPlog 2 P 

(49) 

If only one FIFO is to be used per channel, its capacity must be of H(P-1) 
positions, in order for the total required memory to be HP(P-1). For example, if 
20 H-16 and P=8, according to the preferred method which structure is shown in 
Figure 6, 384 positions of distributed memory would be required, while in the 
case of using on sole FIFO per channel 896 positions would be required. 

FIFO (49) stores temporarily the results obtained by the current sub-module 
25 until the adder (42) depicted in Figure (10) can vacate a position when FIFO 
(52) has a corresponding datum, in which case it will deliver said result to the 
next sub-module. 

Brief description of the contents of the Figures 
30 Figure 1 shows the geometry used as an example with an array composed by 
N numbered elements 1, 2, 3, ...k,...N-1 and N over axis x. The foci are located 
at AR intervals from a minimum distance R 0 . The segment between a focus 
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located in the polar coordinates (R,Q) and the k element of the array located in 
the Cartesian coordinates (x k , 0) which length is L k (R,9) is shown. 

Figure 2 shows the general structure of a generic digital beamforming system 
5 in reception. It entails an array (10) with N elements (11) and a signal 
processing channel per each element. Each channel contains amplification 
circuits and signal conditioning AAS (12) and an A/D digitalizer (13). The 
outputs of the N digitalizers are processed by the beamformer CONF (14). 

10 Figure 3 shows the structure frequently used to execute the coherent 
composition of the signals according to known methods. It is composed by a 
gross delay system (22) and a fine delay system (23) which operation is 
managed from a local controller (24) that uses the contents of a dynamic 
focusing memory (25). The N delayed outputs f h fz-.JN are added in 20 to 

1 5 obtain the resulting r sequence. 

Figure 4 shows the structure of a digital processing channel, preferably used 
for the system that is the subject of this patent. It contains the amplifying and 
signal conditioning element AAS (12), the A/D digitalizer (13), a FIFO memory 
20 (43) that stores temporarily the acquired samples, a sampling clock generator 
(70) that is a fundamental element of the present invention and an enabling 
module of the dynamic focusing and control HFOC (80) . 

Figure 5 shows the initial structure needed to effect the sampling clock 
25 generator (70) with the focaLcorrection coded Q expressed with 1 bit in the 
MEM memory (75). The CNT counter (76) establishes the address to access 
this memory. The sampling clock CK k of channel k is generated by logic 
distributed in the REG-A register (71a), the programmable length shift SHR-A 
(72), the FFA flip-flop (73), the multiplexer MUX (74) and the FFB flip-flop (77). 

30 

Figure 6 shows the initial structure to do the sampling clock generator (70) with 
focal correction codes Q expressed with multiple bits in the MEM memory (75). 
The counter CNT (76) establishes the address to access this memory. The 
sampling clock CKk of channel k is generated by logic distributed in the REG-A 
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register (71a), the programmable length shift SHR-A (72), the programmable 
length shift SHR-B (78) and the set of inverters (75a). 

Figure 7 shows the initial structure to do the sampling clock generator (70) with 
5 focal correction codes Q expressed with one or more bits in the MEM memory 
(75). The counter CNT (76) establishes the address to access this memory. In 
this case a fixed or variable sample number can be inserted between foci, 
according to the contents of the register-counter RCM (71b) controlling the ng 
variable that establishes the length of the shift register (78). Variable ng, of 1 bit 

10 is provided by the CTRL control (90) shown in greater detail in the lower section 
of the figure. This unit is composed by a programmable counter module CMOD 
(91) which contents ant the value of the focal correction code Q form the 
address to one local memory MDM (92). Output g of this memory is inverted by 
(93) to obtain the ng signal. In this case the option to provide the focal offset 

15 code J from the same MEM memory of focal corrections code (75) does exist, 
as shown in a point line or, preferably, from the exterior and through a FIFO 
memory (79). As for the rest, the sampling clock CK k of channel k is produced 
by the logic composed by the REG-A (71a) register, the length shift register 
SHR-A (72) and the length shift register SHR-B (78). The control unit CTRL 

20 provides, also, the cef output used to enable the counting of foci in the HFOC 
(80) unit shown in Figures 4 and 8. 

Figure 8 shows schematically the principle on which the HOFC unit (80) is 
based. This unit does the general control function of the acquisition. It is 

25 composed of a counter CNT-A (82) that is loaded by means of ida with the 
value N A i, it is enabled with the cea signal and it is actuated with the low 
frequency clock cb. The terminal tea count acts on the S input of the FF flip-flop 
(83). Its aj output is delays an amount of master clock c cycles expressed by 
N A2 , by the shift register SHRG (84). The gate AND (81) produces a high level 

30 at the HFD output when its both inputs are on high. The flip-flop (83) is reset to 
0 when its R input is raise when the signal tcf is activated. 

The tcf signal, is the terminal count of the CNT-F counter (86) that can be 
initialized with a specific number of NF foci when the idf signal is activated, its 
35 enabled with the cef signal and actuated with the low frequency cb clock. The 
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CNT-S counter (88) counts the acquired samples, activating its FIN output when 
the N s number programmed is reached. This programming is done by activating 
the ids signal, the counter is enabled with ces and is actuated with the low 
frequency cb clock. 

5 

Figure 9 shows schematically the principle on which the APD unit (60) that 
does the apodization and dynamic aperture functions is based. It is composed 
of a multiplier (61) that obtains the product of input ek times the Ak coefficient 
and delivers the result hk to a multiplexer MUX (65). It also has a CNZ counter 

10 actuated by a cs clock and enabled by the HFD signal, which initial content is 
/Wz, that is, the number of null samples to be provided by this channel to 
implement the dynamic aperture function. When the HFD signal is raised signal 
z is lifted to the output of the FF flip-flop (64), and the output '0' delivered by unit 
(63) by the multiplexer MUX (65) that provides it as output f k is selected. When 

15 the terminal count tc is reached, input R of the FF flip-flop (64) is activated and 
its output z is lowered to a low level prompting the multiplexer MUX (65) to 
select the h k output of the multiplier MULT (61) and is delivered through output 
fk. 

20 Figure 10 shows schematically the modular architecture of the system that is 
the object of this patent. A module (50) contains multiple sub-modules (40), 
each one with various elementary processors PE (30) which input are the 
signals received by the array (10) elements (11). The output of the elementary 
processors is added in (41) and the output of this adder device is used as one 

25 of the terms of the sum in a chain of adders (42) and FIFOs (52). The first sub- 
module substitutes this FIFO by a null input (51). The results provided by the 
adder (42) of the last sub-module are temporarily stored in a FIFO (56) to be, 
then added themselves in (58) to the results of a preceding module that have 
been stored in FIFO (57) in adder (58). Additionally the module contains a 

30 memory MJ (59) that provides the focal offset codes J to the different 
elementary processors PE (30). 

Figure 11 shows schematically, the principle on which the preferred 
embodiment of the adder (41) is based, with an example for 8 channels where 
35 multiple FIFOs (43) store temporarily the samples obtained at each elementary 
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processor ek. The FIFOs (43) outputs are added by pairs in the adders (44), 
and the result stored in the FIFOs (47). Similarly, the outputs of said FIFOs are 
added in (48), and the final results are stored in FIFO (49) which output r 
represents the sum of the samples e 1 to e 8 delivered by the 8 elementary 
5 processors that are a sub-module in this example. 

Figure 12 shows a diagram of an example of a module (200) of the system that 
is the object of the present patent. It has two connectors (120a) and (120) b 
through which the analogic signals to be distributed by the buses (121a) and 

10 (121b) are entered to the double A/D converters (101). The samples provided 
by these converters are introduced by (106) in a FPGA (102) that does the 
digital processing functions corresponding to a sub-module (100), amongst 
which are the generation of the sampling clocks (105) and the provision of the 
results to the next sub-module via the bus (103). The last sub-module in the 

15 chain provides the results obtained to a FPGA (111) by the bus (104) to be 
combined with the preceding results of another module by bus (113), the 
subsequent results provided by buses (112) or (114). These buses are present 
in the connector (110). Also the FPGA (111) delivers control signals (115a) and 
(115b) to the amplifier and analogic conditioning circuits via connectors (130a) 

20 and (130b) respectively. 

Figure 13 shows the experimental measures done on the module described in 
Figure (12), showing the absolute error in ns for each of the 32 channels. 

25 

EXAMPLE OF AN EMBODIMENT OF THE INVENTION 

In all the cases described the sampling clock for each channel is obtained using 
very few hardware resources, essentially registers, counters and a memory 
used to store the focal correction codes, and the signal acquisition is achieved 
30 with a high temporal resolution. 

Due to the high efficiency of the methodology described in what pertains to the 
use of hardware resources, and contrary to other embodiments, the PFC and 
the VPFC allow for the integration of multiple processing channels in one sole 
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VLSI circuit, with other functions as described, even in devices of general use 
such as the FPGAs. 

In a demonstration embodiment of this invention, a module contains 4 sub- 
5 modules, each of them with 8 elementary processors. This small system can 
operate with an active aperture of up to 32 elements in linear or sectorial scan 
with disperse aperture and arrays of 128 elements. A system that operates with 
dense or active apertures of up to 128 elements is configured with 4 modules. 
The system is scalable, adding 32 active elements for each new module 
10 installed. 

Figure 12 shows graphically the module's configuration done in a printed circuit 
card of 220x100 mm. the circuits corresponding to the analogic pre-processing 
of the signals with a possible excitation of the transducers, as well as the post- 
15 processing circuits, that in this particular embodiment are based on an 
integrated computer to which the Interface and Control Unit (ICU) is connected 
to the system described have been omitted. 

Module (200) receives the analogic signals from 32 channels provided by the 
20 amplifier and conditioning circuits via the connectors (120a) and (120b). It 
contains four sub-modules (100) and one combiner done in the FPGA (111). 
The combiner receives signals processed by the bus (113) in a preceding 
module in the chain and sends the results to the following one by means of bus 
(114). The parameters are programmed from the ICU via the control bus (112), 
25 via which also the last module in the chain delivers the final results. Also, the 
FPGA (111) delivers a set of digital signals via the connectors (130a) and 
(130b) that are used by the analogic processing unit. This structure faithfully 
reflects the architecture shown in Figure 10. 

30 The analogic signals introduced by the connectors (120a) and (120b) are 
distributed to the sub-modules via buses (121a) and 121b respectively, where 
they are digitalized by the A/D converters (101). The digital outputs of the A/D 
converters are introduced in the FPGAs (10) via the buses (106). The sampling 
clock of said converters (105) is generated in the FPGAs (102) with the 

35 Progressive Focal Correction technique described in this specification. 
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The sub-modules are connected amongst themselves via the bidirectional 
buses (103). Together with bus (104), that connects the last sub-module to the 
combiner (111), these buses are used in one direction for the programming of 
5 parameters from the control one (112) and in the opposite direction to provide 
results to the sub-modules that are next in the chain or to the combiner (111). 

The combiner does the coherent summation of the results obtained at the 
current module via the bus (104) to those results received from a preceding 
10 module by bus (113), providing the results of this process via bus (114) to the 
next module, or via bus (112) to the Interface and Control Unit. Either option is 
programmable. 

In this embodiment the following components have been used: 16-10 bits each- 
15 , type AD9218 double converters of Analog Devices have been used for the 
devices (101), 4 type XC2S200 FPGAs from Xilinx Inc. for devices (102) and 
one XC2S150 from Xilinx Inc., for device (111). 

The system clock operates at 40 MHz, a speed at which the communication 
20 between sub-modules, modules and the ICU operate synchronically. Internally, 
part of the circuits also operate at this frequency or at a fraction that is 
determine by the sampling frequency, however, in the FPGAs (102), it is 
multiplied by 4 times the frequency of the system clock to obtain the 160 MHz 
master clock. The temporal resolution is then Tx = 6.25 ns and the sampling 
25 error is theoretically delimited to Id < 3.125 ns. This allows the processing of 
analogic signals of up to 10 MHz with p £ 6 and up to 20 MHz with p> 8, where 
p is the relationship between the signal's period and the temporary resolution 
that affects the dynamic range and contrast of the images. 

30 Each FPGA (102) contains the necessary circuits to implement all the digital 
processing of 8 channels, including: 

■ Apodization, with 0 < A k < 255/256 coefficients or disabled (A k = 1, Vk) 

■ Dynamic aperture, where 0 < MZ k < 2 12 

35 ■ Beam deflection with a fixed focus, with 0 < NA k < 2 14 
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■ Dynamic focusing according to the Progressive Focal Correction theory 
as described above, with 1 bit 4096 focal correction codes per each 
channel 

■ Programming the number of samples to be acquired, where 2 <Ns< 2 17 
5 ■ Programming the number of samples between foci, where 4<m<64 

■ Programming of the number of foci per acquisition, where 0<NF< 2 12 

■ Programming of the output sampling frequency, where 2.5<fs<40 MHz 

■ Circuits managed for energy efficiency and enabling of functions 

■ Parameter programming and processing of results circuits 

10 ■ Circuits to verify the module, integrating a digital signal generator 

The combiner contained in the FPGA (111) has circuits to do the following 
tasks: 

15 ■ Module management, transferring parameter values to the sub-modules 
via (104) and (103) 

■ Compilation of results obtained from the preceding module via (113) 

■ Composition of said results with those of the module, received via (104) 

■ Provision of the combined results to the following module via (114) 
20 ■ Provision of the final results to the ICU via bus (112) 

■ Circuits to compensate signal offset 

■ Circuits to control the analogic pre-processing unit 

Figure (13) shows experimental measures relative to the absolute error in ns (in 
25 the ordinates) of the sampling instant in each of the 32 channels (in the 
abscissas) for a focus located at a 50 mm distance, with a 15° deflection angle, 
propagation speed of c = 5900 ms' 1 , and a 5 MHz linear array with elements 
uniformly distributed at intervals of d = 0.8 mm. 

30 The maximum theoretical error is of 3.125 ns, which is verified in 28 of the 32 
channels. The greatest error in the 4 channels that exceeds the theoretical 
value is due to differences in the track layouts in the printed circuit and in the 
FPGAs of this experimental model, and is not a significant error at any rate. The 
average square error, which is more determining of the quality of the image, is 

35 of 1.9 ns. 
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Evidently, improvements may be introduced to the embodiment described 
herein, since there are devices that have a greater integration capacity. Also, 
other modifications that affect the type of A/D conversion, the distribution of 
5 devices, etc., can be introduced, since the model described is only submitted as 
a demonstrative example, without said modifications entailing a substantial 
change from the tenets of the methodology and system described iin the 
present specification. 

10 To apply the method, a case of sectorial ultrasonic image in pulse-echo with a 
linear transducer of N elements centered at origin and size D aperture can be 
considered. The aim is to cover a -45° £ 6 £ 45° sector with the image. 

If f R is the frequency of the received signal, in order to verify with sufficient 
15 margin the Nyquist criterion, a nominal sampling frequency of f s = 4 f R is 
selected. Also, to maintain a high dynamic range and contrast in the image p = 
fx/ f R = 16 is selected, where f x is the master clock frequency. 

Considering the round trip travel path of the ultrasonic pulse, the interval 
20 between samples is: 

(50) &R=cT s /2 

25 

Substituting the AR value given by Equation 1, when m = 1 the result obtained 

is: 




(51) 

30 When m samples are inserted between foci, the value of v = 4 m. Applying now 
Equation 18, with lx k l ^ D/2, where D is the size of the aperture, the minimum 
distance to apply the Progressive Focal Correction technique with complete 
aperture is: 



51 



D 2 



mv -2a 



rCOS0 



amv -a 



(52) 



5 

Where a =2 b - 1, where b is the number of bits with which to express the focal 
correction codes, provided that mv > 2a. 

Table II shows the results of Equation (52) for b =1, 2 and 3 based on the 
10 (fixed) number of samples between foci. 

It must be observed that the value Ro/D represents the minimum value of the 
number F# that can be used to implement the operation using the dynamic 
aperture function. Generally, a value of F# mjn ^ 1 that is obtained with most of 
15 the combinations is accepted as valid. 

Table II 

F# min = Rq/D to apply the PFCT with complete aperture (-45°< 0 < 45°) 



M 


b=1 


b=2 


b=3 


1 


0.6 






2 


0.8 


0.4 




4 


1 


0.6 


0.4 


8 


1.4 


0.9 


0.6 



20 
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